Constraining the Quintessence equation of state 
with Snia data and CMB peaks 
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ABSTRACT 

Quintessence has been introduced as an alternative to the cosmological constant scenario 
to account for the current acceleration of the universe. This new dark energy component allows 
values of the equation of state parameter w^ > — 1 , and in principle measurements of cosmological 
distances to Type la supernovae can be used to distinguish between these two types of models. 
Assuming a flat universe, we use the supernovae data and measurements of the position of 
the acoustic peaks in the Cosmic Microwave Background (CMB) spectra to constrain a rather 
general class of Quintessence potentials, including inverse power law models and recently proposed 
Supergravity inspired potentials. In particular we use a likelihood analysis, marginalizing over 
the dark energy density VIq, the physical baryon density Vlbh^ and the scalar spectral index n, 
to constrain the slopes of our Quintessence potential. Considering only the first Doppler peak 
the best fit in our range of models gives vJq ~ —0.8. However, including the SnIa data and 
the three peaks, we find an upper limit on the present value of the equation of state parameter, 
— 1 < Wq < —0.93 at 2(7, a result that appears to rule out a class of recently proposed potentials. 

Subject headings: Cosmic Microwave Background Anisotropy, Cosmology 



1. Introduction 

Observations of distant type la supernovae 
(Perlmutter et al. 1999; Riess et al. 1999) and 
small angular scale anisotropics in the Cosmic Mi- 
crowave Background (CMB) (De Bernardis et al. 
2000; Balbi et al. 2000; Netterfield et al. 2001; 
Pryke et al. 2001) suggest that the universe is 
dominated by a large amount of dark energy with 
a negative equation of state parameter w. One 
obvious explanation would be the presence for all 
time of a cosmological constant with w — —1, al- 
though there is no satisfactory reason known why 
it should be so close to the critical energy den- 
sity (for a general review see Sahni & Starobinsky 
2000). An alternative proposal introduces a new 
type of matter and is called 'Quintessence' (Cald- 
well et al. 1998). Assuming that some unknown 
mechanism cancels the true cosmological constant, 



this dark energy is associated with a light scalar 
field Q evolving in a potential V{Q). The equation 
of state parameter of the Q component is given by 
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and it is a function of time. According to the 
form of V{Q) the present value of wq is in the 
range v/q > —1. The temporal dependence 
of WQ implies that high red-shift observations 
could in principle distinguish between KCDM 
and QCDM models (Maor et al. 2000; Huterer 
& Turner 2000; Alcaniz & Lima 2001; Benabed & 
Bernardeau 2001; Cappi 2001; Weller & Albrecht 
2001). Moreover, a number of authors have re- 
cently pointed out that the position of the CMB 
peaks could provide an efficient way to constrain 
Quintessence models (Kamionkowski & Buchalter 
2000; Croocks et al. 2000; Doran et al. 2000). 

In this paper we use the supernovae sample 
of Perlmutter et al. (1999) and the recent mea- 



surements of the location of the CMB peaks (De 
Bernardis et al. 2001) to determine new limits on 
the Quintessence equation of state. Our study is 
similar in approach to an earlier analysis by Ef- 
stathiou (2000). We consider a general class of 
potentials parametrized in such a way that we can 
control their shape, and apply a likelihood analysis 
to find the confidence regions for the parameters of 
the potential and the best value for the fractional 
Quintessence energy density fig. The constraints 
which emerge are different if we analyze the data 
separately. In particular the position of the first 
Doppler peak prefers a Quintessence model with 
v/q ~ —0.8 for the prior Q.q = 0.7 in agreement 
with Baccigaluppi et al. (2001), while the analy- 
sis of all the CMB peaks and Snia gives an upper 
value for the equation of state, Wg < —0.93 at 2(j 
for these class of models. This limit is stronger 
than those previously obtained (Perlmutter et al. 
1999; Efstathiou 2000; Amendola 2001; Balbi et 
al. 2001), Wg < —0.6 at 2cr, simply because we are 
making use of the new improved CMB data. An 
obvious consequence of this result is that in these 
class of models, for them to succeed the scalar field 
dynamics has to produce effects similar to pure 
vacuum energy and in this case it is unlikely that 
Quintessence can be distinguished from a cosmo- 
logical constant (see also Maor et al. 2000). 

2. Quintessence equation of state 

The scalar field dynamics is described by the 
Klein-Gordon equation 
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where pm and pr arc the matter and radiation en- 
ergy densities respectively. It is well known that 
for a wide class of potentials, Eq. (2) possesses at- 
tractor solutions (Steinhardt et al. 1999). In this 
regime the kinetic energy of the field is subdomi- 
nant allowing wq to become negative. The present 
value of Wq depends on the slope of the poten- 
tial in the region reached by the field. Actually if 
the Quintessence field rolls down a very flat region 



(Barreiro et al. 2000) or if it evolves close to a min- 
imum (Albrecht & Skordis 1999, Brax & Martin 
1999, Copeland et al. 2000) the equation of state 
parameter varies in the range — 1 < Wq < —0.8. 
On the other hand models like the inverse power 
law potential (Ratra & Peebles 1988; Zlatev et 
al. 1999) require larger values of Wq. A general 
potential which can accomodate a large class of 



scenarios is: 
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where k — \/8TrG and M is fixed in such a way 
that today pg = Pc^Q, where pc is the critical en- 
ergy density. For (3 = Eq. (4) becomes an inverse 
power law, while for (3 = 2 we have the SUGRA 
potential proposed by (Brax & Martin 1999). For 
a = 0, (3 = 1 and starting with a large value of 
Q, the Quintessence field evolves in a pure expo- 
nential potential (Ferreira & Joyce 1998). We do 
not consider this case further since it is possible 
to have a dark energy dominated universe, but at 
the expense of fine timing for the initial condi- 
tions of the scalar field. Larger values of /3 mimic 
the model studied by Copeland et al. (2000). For 
a, (3 ^ the potential has a minimum, the dynam- 
ics can be summarized as the following. For small 
values of /3 and for a large range of initial condi- 
tions, the field does not reach the minimum by the 
present time and hence Wq > —1. For example, 
if the Quintessence energy density initially domi- 
nates over the radiation, the Q field quickly rolls 
down the inverse power law part of the potential 
eventually resting in the minimum with wq ^ — 1 
after a series of damped oscillations (Riazuelo & 
Uzan 2000). This behaviour however requires fine 
tuning the initial value of Q to be small. On the 
other hand, this can be avoided if we consider large 
values of a and /3 (Fig. la). In these models the 
fractional energy density of the Quintessence field, 
f2g , is always negligible during both radiation and 
matter dominated eras. In fact, for small initial 
values of Q, V(Q) acts like an inverse power law 
potential, hence as Q enters the scaling regime its 
energy density is subdominant compared to that of 
the background component. Therefore nucleosyn- 
thesis constraints (Bean et al. 2001) are always 
satisfied and there are no physical effects on the 
evolution of the density perturbations. The main 



consequence is that for a different value of Wq the 
Universe starts to accelerate at a different red-shift 
(Fig. lb). This implies that different values of a 
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Fig. 1. — In (a) the evolution of wq against the 
red-shift is plotted for different values of a and 
p. In (b) the behaviour of the deceleration pa- 
rameter, q, is plotted against the red-shift. The 
acceleration starts {q < 0) earlier for models with 
an equation of state close to that of a true cosmo- 
logical constant. 

and (3 lead to a different luminosity distance and 
angular diameter distance. Consequently by mak- 
ing use of the observed distances we may in prin- 
ciple determine an upper limit on Wq, potentially 
constraining the allowed shape of the Quintessence 
potential (Huterer & Turner 2000). 

3. CMB peaks 

The CMB power spectrum provides informa- 
tion on combinations of the fundamental cosmo- 
logical quantities. The position of the Doppler 
peaks depends on the geometry of the Universe 
through the angular diameter distance, although 
the amplitude of the peaks are sensitive to many 
different parameters. The important point for us 
is that in general the Quintessence field can con- 
tribute to the shape of the spectrum through both 
the early integrated Sachs- Wolfe effect (ISW) and 
the late one (Hu et al. 1997). The former is 
important if the dark energy contribution at the 
last scattering surface (LSS) is not neghgible (Sko- 
rdis & Albrecht 2000, Barreiro et al. 2000) or in 
non-minimally coupled models (Amendola 2000, 



Perrotta et al. 2000, Baccigalupi et al. 2000), 
whereas the late ISW is the only effect in models 
with f^Q ~ at LSS (Brax et al. 2000). However, 
as has recently been demonstrated an accurate de- 
termination of the position of the Doppler peaks 
is more sensitive to the actual amount of dark en- 
ergy (Doran et al 2000). To be more precise, the 
multipole of the ?7T,-th peak is Im = rnlsh , where Ish 
is proportional to the angular scale of the sound 
horizon at LSS. In a flat universe Ish is given by: 
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where Cs is the mean sound velocity and tq, tis are 
the conformal time today and at last scattering re- 
spectively. However, physical effects before recom- 
bination can shift the scale of the sound horizon at 
different multipoles, resulting in a better estimate 
for the peak positions being given by: 



Im = Ishim - SI - 6l,n), 



(6) 



where SI is an overall shift (W. Hu et al. 2000) and 
Sim is the shift of the m-th peak. These correc- 
tions depend on the amount of baryons flbh^, on 
the fractional quintessence energy density at last 
scattering {^q) and today {^q), as well as on the 
scalar spectral index n. Recently, analytic formu- 
lae, valid over a large range of the cosmological 
parameters, have been provided to good accuracy 
for SI and Sim (Doran & Lilley 2001). Of crucial 
importance is the observation that the position of 
the third peak appears to remain insensitive to 
other cosmological quantities, hence we can make 
use of this fact to test dark energy models (Doran 
et al. 2001). 

4. Likelihood analysis and results 
4.1. Constraints from supernovae 

We want to constrain the set of parameters a, 
P and Oq confined in the range: a G (1,10), 
/3 e (0,10) and flq e (0,1), subject to the as- 
sumption of a flat universe. We use the Snia data 
fit C of Perlmutter et al. (1999), that excludes 4 
high redshift data points. The magnitude-redshift 
relation is given by: 



m(z) = M + 5logVLiz,a, P^Qq), 



(7) 



where M is the 'Hubble constant free' absolute 
magnitude and Vl = HodL{z) is the free- Hubble 
constant luminosity distance. In a flat universe 
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where tq is the conformal time today and t{z) is 
the conformal time at the red-shift z of the ob- 
served supernova. Both of these quantities are 
calculated solving numerically Eq. (2) and Eq. (3) 
for each value of a,/3 and fig. In Ai we neglect 
the dependence on a fifth paramenter [a in Perl- 
mutter et al. 1999) and assume it to be 0.6, the 
Perlmutter et al. (1999) best value. We then ob- 
tain a gaussian likelihood function £'^"(a, (3, ^q), 
by marginalizing over M. In Fig. 2a we present the 
one-dimensional likelihood function normalized to 
its maximum value for Qq. There is a maximum 
at rig = 1, in agreement with Efstathiou (2000). 
In Fig. 3a we present the likelihood contours in the 
a — p parameter space, obtained after marginal- 
izing over i^Q. Note that all values are allowed 
at the 2ct level. The confidence regions for the 
Snia data correspond to Quintessence models with 
Wq < —0.4 for Hq — 0.6, an upper limit that 
agrees with both Perlmutter et al. (1999) and Ef- 
stathiou (2000). 
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Fig. 2. — Fractional Quintessence energy density 
likelihoods, (a) for SnIa, (b) for the combined 
CMB peaks and (c) for the combined data sets. 
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Fig. 3. — Likelihood contour plots for SnIa, I, II 
and III acoustic peaks. The blue region is the 68% 
confidence region while the 90% is the light blue 
one. For the SnIa the white region correspond 
to 2a. The position of the third CMB acoustic 
peak strongly constrains the acceptable parameter 
space. 



4.2. Constraints from Doppler peaks and 
SnIa 

We now compute the position of the three 
Doppler peaks li,l2 and I3 using Eq. (6). In ad- 
dition to the parameter space used in the super- 
novae analysis we consider the physical baryon 
density and the scalar spectral index varying re- 
spectively in the range Qth^ S (0.018,0.026) and 
n e (0.9,1.1). The Hubble constant is set to 
h = 0.70 in agreement with the recent HST ob- 
servations (Freedman et al. 2000). The pre- 
dicted peak multipoles in the CMB are then com- 
pared with those mesured in the BOOMERANG 
and DASI spectra (De Bernardis et al. 2001). 
Note, that the third peak has been detected in 
the BOOMERANG data but not in the DASI 
data. Furthermore the authors of De Bernardis 
et al. (2001), with a model independent anal- 
ysis, estimated the position of the peaks accu- 
rately at la. However because the errors asso- 
ciated with the data are non Gaussian, to be con- 
servative we take our la errors on the data to 



be larger than those reported in De Bernardis et 
al. (2001), so that our analysis is significant up 
to 2a. We then evaluate a gaussian likelihood 
function C^^"'^'^ {a, j3,D.Q,Vl},h? ,n). The combined 
one-dimensional likelihood function for the peaks 
is shown in Fig. 2b, where we find VIq = 0.69± q-^q . 
The likelihood for all the data sets combined is 
shown in Fig. 2c, where we find fig = 0.75 ± [J^g. 
These results are in agreement with the analysis 
of Efstathiou (2000), Netterfield et al. (2001) and 
Baccigalupi et al. (2001). 
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Fig. 4. — One-dimensional likelihood for n and 

The likelihood functions, combining all the data 
for the CMB peaks, for the scalar spectral in- 
dex and the physical baryon density are shown in 
Fig. 4. Since the dependence of the peak multipoles 
on r2{,/i^ and n is small, it is not possible to obtain 
some significant constraints on these cosmologi- 
cal parameters using the location of the Doppler 
peaks. In Fig.3b-3d we plot the two-dimensional 
likelihood function in the plane a — [3 for each 
peak, obtained after having marginalized over Jig, 
Vlbh? and n. Their shape reflects the accuracy in 
the estimation of the position of the peaks. Ac- 
tually the first one is very well resolved, while we 
are less confident with the location of the second 
and third peak. Therefore their likelihoods are 
more spread and fiat in the a — (3 plane. The la 
confidence contour (Fig. 3b) for the first acoustic 
peak constrains the slopes of our potential in the 
range: 3 < a < 10 and 1 < /3 < 3. In par- 
ticular the likelihood has a maximum at a = 9 




Fig. 5. — Two-dimensional likelihood for Snia and 
CMB with 1 (dark blue) and 2cr (light blue) con- 
tours. 



and [3 = 2, corresponding to an equation of state 
w% = —0.8 for 51q = 0.7, in agreement with the 
recent analysis by Baccigalupi et al. (2001). How- 
ever, the second and third peaks constrain a region 
where the equation of state is compatible with the 
cosmological constant value. Therefore the effect 
of including all the data in the likelihood anal- 
ysis is to move the constraint from models with 
Wq ~ —0.8 to models with an equation of state 
Wg '--^ — 1. As we can see in Fig. 6 the values of 
a and /3, allowed by the likelihood including all 
the data (Fig. 5), correspond to our models with 
values of Wg in the range — 1 < w% < 0.93 at 
2a for our prior probability f2g = 0.75. The rea- 
son for such a strong constraint is due to the as- 
sumed accurate determination of the third peak, in 
that it is insensitive to pre-recombination effects. 
In particular peak multipoles are shifted toward 
larger values as Wq approaches the cosmological 
constant value. This is because, in models with 
Wq r^ —1 the universe starts to accelerate earlier 
than in those with Wq > — 1, consequently the dis- 
tance to the last scattering surface is further and 
hence the sound horizon at the decoupling is pro- 
jected onto smaller angular scales. Since the loca- 
tion of the third peak inferred by De Bernardis et 
al. (2001) is at I3 = 845± ^^ , values of w^ 1 fit 
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Fig. 6. — Equation of state parameter against a 
and /3 for l^g = 0.75. The 2a contours correspond 
to models with vlk '^^ — 1- 



this multipole better than models with Wg > — 1. 
However we want to point out that at \a the po- 
sition of the first peak is inconsistent with the po- 
sition of the other two. A possible explanation of 
this discrepancy is that the multipoles I2 and Iz 
are less sensitive to small shift induced by the de- 
pendence on fibft,^ and n. Therefore we can obtain 
a different constraint on the dark energy equation 
of state if we consider the peaks individually. 

5. Conclusions 

The location of the sound horizon is very sen- 
sitive to the dark energy contribution. Due to 
the strong degeneracy in the shape of the CMB 
spectrum, a certain class of Quintessence models 
can be better constrained using only the acoustic 
peaks. We have applied a likelihood analysis to 
constrain the shape of the Quintessence potential, 
based on both the supernovae type la data and 
the positions of the CMB peaks. Assuming a flat 
space-time and making use only of the position 
of the first Doppler peak we find the best fit for 
models w% ~ —0.8 for fig = 0.7 prior value. The 
combined analysis, including all three peaks and 
Snia, gives the best fit for ?Lq = 0.75 ± °°g. We 
have found in particular that the determination 
of third peak in the BOOMERANG data hmits 
the equation of state parameter at 2a in the range 
— 1 < Wq < —0.93 for rig with this prior value. 



This has an important implication for minimally 
coupled Quintessence models. Actually they must 
behave similarly to a cosmological constant, there- 
fore inverse power law is disfavoured. In fact, an 
equation of state parameter iiPq ~ — 1 implies the 
Quintessence field is undergoing small damped os- 
cillations around a minimum or evolving in a very 
flat region of the potential. For these reasons mod- 
els like the double exponential potential (Barrciro 
et al. 2000) or the single modified exponential 
potential (Skordis & Albrecht 2000) pass this con- 
straint, even though they are not included in our 
analysis. Another important caveat is that this 
study does not take into account Quintessence sce- 
narios where the contribution of the dark energy 
density in radiation or early in matter dominated 
eras is not negligible. In such a case we would 
have to take into account physical effects not only 
on distance measurements, but also on the struc- 
ture formation process itself. These models and 
the non-minimally coupled ones therefore could 
yet be distinguished from a pure KCDM model. 
We still require a more complete analysis to un- 
derstand the nature of the dark energy, but this 
paper points out that it is possible to constrain 
certain classes of models far more than was previ- 
ously reafised. 

We are very grateful to the referee for mak- 
ing a number of constructive comments and to L. 
Amcndola, A. Liddle, N.J. Nunes and J. Weller 
for useful discussions and comments. 
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